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Introduction 

The  specific  purpose  of  this  research  contract  is  to  study  the  scattering  of 
focused  electromagnetic  beams  by  lossy  targets.  This  objective  has  been  duly  accom¬ 
plished,  and  as  a  result  of  an  elaborate  study  of  the  subject  in  a  Ph.D  dissertation.  The 
broader  objective  of  this  research  study  is  the  general  application  of  time  domain  com¬ 
putational  techniques  in  electromagnetic  theory.  That  general  objective  is  also  accom¬ 
plished  in  numerous  papers  published  and  presented  in  the  area  of  time  domain 
methods. 

Accomplishments 

The  general  objective  of  this  contract  is  the  study  of  the  application  of  time 
domain  techniques  in  electromagnetic  computation.  We  have  used  the  specific  goal  of 
understanding  the  interaction  of  focused  beam  with  lossy  media  to  isolate  many  impor¬ 
tant  issues  associated  with  time  domain  computations  yet  can  also  be  applied  to  other 
areas,  notably,  microstrips  and  MMIC  related  computations.  Those  issues  are: 

(1)  Absorbing  Boundary  Conditions 

In  the  computation  of  open  systems,  the  finite  difference  or  finite  element 
meshes  must  be  terminated  by  an  artificial  absorbing  boundary  condition  which 
simulates  a  non-reflective  boundary.  Since  this  is  such  an  important  issue  with 
time  domain  computation,  the  absorbing  boundary  is  indeed  a  part  of  differential 
equation  approach  to  scattering  problems.  There  are  many  ways  an  absorption 
condition  can  be  implemented.  Some  absorbing  conditions  are  easy  to  apply  but 
do  not  absorb  well  and  some  absorbing  conditions  absorb  well,  but  are  not  so 
easy  to  apply. 

In  our  investigation,  we  have  invented  a  procedure  termed  “Superabsorp¬ 
tion.”  This  method  can  improve  any  existing  absorbing  boundary  condition  by 
an  order  of  magnitude  with  very  simple  programming  procedures.  It  is  indeed  a 
breakthrough  in  the  study  of  absorbing  boundary  conditions. 

(2)  Propagation  Along  Interfaces 

The  time  domain  investigation  of  wave  propagation  along  dielectric  inter¬ 
faces  is  another  important  issue  in  TD  computation.  A  typical  such  problem  is  a 
microstrip  at  high  frequencies  when  the  quasistatic  results  fail.  We  have  used  the 
time  domain  method  to  investigate  the  performance  of  microstrips  and  other 
MMIC  components. 

(3)  Penetration  of  Focused  Beam  into  Lossy  Media 

The  penetration  of  focused  beam  is  an  effort  to  find  out  whether  the  tem¬ 
perature  of  a  body  can  be  raised  in  the  interior  without  overheating  the  exterior. 
This  is  an  issue  most  interested  by  those  interested  in  biological  effects  of 
microwaves  and  hyperthermia.  Our  computation  shows  that  it  is  possible  to 
make  a  hot  spot  using  focused  beam. 


(4)  Computations  Involving  Superconductive  Material 

There  are  many  experiments  now  going  on  regarding  the  application  of 
superconductors  in  microwave  and  mm-waves.  Indications  are,  one  day  this 
material  will  be  used  in  many  devices.  Any  serious  application  of  this  material 
will  involve  computations.  How  should  we  treat  this  material?  Is  a  superconduc¬ 
tor  the  same  as  a  perfect  conductor?  These  questions  are  answered  in  our  paper 
on  “Electromagnetics  of  Superconductors.” 

In  the  following  pages,  we  have  listed  the  abstract  of  the  papers  that  have  been 
published  or  presented  under  the  sponsorship  of  this  contract. 
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TJnimoment  Method  for  Electromagnetic  Wave  Scattering 

K.  K.  Mei 

Department  of  Electrical  Engineering 
University  of  California 
Berkeley,  CA  94720,  USA 


Abstract-A  simple  problem  in  2-dimensional  space  is  used  as  an  example  for  the  anolira 
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oadmg  of  equations,  sparse  matrix  inversion  and  boundary  matching.  The  objective  of 
this  paper  is  to  help  uninitiated  readers  understand  and  apply  this  method!  J 
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Tutorial _ 

Point-Matched  Time  Domain  Finite  Element 
Methods  for  Electromagnetic  Radiation 

and  Scattering 

ANDREAS  C.  CANGELLARIS,  member,  ieee,  CHUNG-CHI  LIN,  member,  ieee,  and  KENNETH  K.  MEI,  fellow,  ieee 


Abstract — Direct  time  domain  computation  of  Maxwell’s  differential 
equations  will  soon  become  a  practical  technique  because  of  the 
availability  of  supercomputers.  Tbe  principal  methods  used  in  time 
domain  computations  and  the  supporting  theories  are  presented.  The 
point-matched  finite  element  method  is  chosen  as  the  main  feature  of  this 
presentation,  which  includes  the  discretization  of  equations,  conforming 
mesh  generation,  dielectric  and  metallic  interfaces,  numerical  stability 
and  simulation  of  radiation  conditions.  Numerical  results  of  scattering  of 
Gaussian  pulses  are  presented  in  time  sequences. 


SS»o«S0N  CHARACTERIST1CS  of 

niC*R<iSTR,P  UNES  ,N  THE  vicinity 
OF  A  COPLANAR  GROUND 


Indexing  terns:  Microwave  circuits  and  systems.  Microstrip 

gr0Und  on  thed«fW«on  charac- 
!fnm  /  ^Unn8  nucrostnP  is  Studied  using  the  time- 

JJSJS’S  erUX  meth0i  ChanS“  in  characteristic 
unpedana  and  propagation  constant  from  those  of  ideal 
nucrosmp  are  obtained.  Further,  an  imaginary  pan  of  the 

'mpe*aace  and  a  non-negligible  attenuation 
constant  are  detected  at  high  frequencies. 


IEEE  TRANSACTIONS  ON  MICROWAVE  THEORY  AND  TECHNIQUES,  VOL  36.  NO.  2,  FEBRUARY  1988 


263 


Calculations  of  the  Dispersive  Characteristics 
of  Microstrips  by  the  Time-Domain 
Finite  Difference  Method 

XIAOLEI  ZHANG,  JIAYUAN  FANG,  KENNETH  K.  MEI,  fellow,  ieee,  and  YAOWU  LIU 


Abstract — The  dispersive  characteristics  of  microstrips  have  been  in¬ 
vestigated  by  many  authors  (1)— [4|  using  various  numerical  and  empirical 
methods,  'lliose  results  showed  a  lack  of  agreement  with  each  other,  and 
the  true  dispersive  characteristics  of  microstrips  still  need  to  be  identified. 
In  this  paper,  a  direct  time-domain  finite  difference  method  is  used  to 
recharacterize  the  microstrip.  Maxwell's  equations  are  discretized  both  in 
time  and  space  and  a  Gaussian  pulse  is  used  to  excite  the  microstrip.  The 
frequency-domain  design  data  are  obtained  from  the  Fourier  transform  of 
the  calculated  time-domain  field  values.  Since  this  method  is  completely 
independent  of  all  the  above-mentioned  investigations,  the  new  results  can 
be  considered  as  an  impartial  verification  of  the  published  results. 

The  comparison  of  the  time-domain  results  and  those  from  the 
frequency-domain  methods  has  shown  the  integrity  of  the  time-domain 
computations.  Since  this  method  is  very  general  and  can  be  applied  to 
model  many  other  microwave  components,  its  success  in  the  microstrip 
problem  is  an  important  step  toward  its  general  application. 
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Time-Domain  Finite  Difference  Approach  to 
the  Calculation  of  the  Frequency-Dependent 
Characteristics  of  Microstrip  Discontinuities 

XIAOLEI  ZHANG  and  KENNETH  K.  MEI,  fellow,  ieee 


Abstract  — The  frequency -dependent  characteristics  of  the  microstrip 
discontinuities  have  previously  been  analyzed  using  several  full-wave  ap¬ 
proaches.  The  time-domain  finite  different  (TD-FD)  method  presented  in 
this  paper  is  another  independent  approach  and  is  relatively  new  in  its 
application  for  obtaining  the  frequency-domain  results  for  microwave 
components  [26].  The  purpose  of  this  paper  is  to  establish  the  validity  of 
the  TD-FD  method  in  modeling  circuit  components  for  MMfC  CAD 
applications. 
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PENETRATION  OF  A  FOCUSED  ELECTROMAGNETIC  PULSE 
INTO  A  LOSSY  MEDIUM 

H.  C.  Chang  and  K.  K.  Mei 

Department  of  Electrical  Engineering  and  Computer  Sciences 
and  the  Electronics  Research  Laboratory, 

University  of  California,  Berkeley,  CA  94720 

Interaction  of  electromagnetic  beam  with  a  lossy  interface  is  studied  using  a 
time  domain  finite  element  method.  A  plane  wave  of  Gaussian  pulse  is  transmitted  towards  a 
2-D  parabolic  reflector,  which  causes  a  high  intensity  field  around  the  focus,  as  the  rest  of  the 
wave  passes.  When  this  focusing  field  is  directed  towards  an  air-lossing  medium  interface,  the 
time  sequence  of  the  penetration  of  the  fields  into  the  lossy  medium  can  be  observed.  If  the 
heat  conductivity  and  specific  heat  of  the  medium  is  known,  the  temperature  profile  of  the 
lossy  medium  near  the  area  of  penetration  can  also  be  calculated  in  time  sequence.  This  talk 
presents  the  basic  numerical  techniques  and  some  recently  obtained  results  of  calculation. 
Further  development  of  this  project  can  have  interesting  biomedical  applications. 
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B4-7  TO  BE  OR  NOT  TO  BE  IN  THE  TIME  DOMAIN 

1620  Kenneth  K.  Mei 

Department  of  Electrical  Engineering  and  Computer  Sciences 
University  of  California 
Berkeley,  CA  94720 

Maxwell's  differential  equations  in  the  frequency  domain  are  ellip¬ 
tical,  which  results  in  classical  boundary  value  problems.  In  the  time 
domain  they  are  hyperbolic,  the  solutions  of  which  march  in  time.  Nu¬ 
merically,  the  frequency  domain  differencial  equations  require  solutions 
of  matrices;  the  time  domain  differential  equations,  if  explicit,  are 
spared  from  such  endeavor.  Before  the  coming  of  super  computers,  almost 
all  electromagnetic  problems  were  solved  in  the  frequency  domain.  Now, 
there  is  a  choice.  This  calk  enumerates  the  pros  and  cons  of  computations 
in  either  domain.  Depending  on  the  objectives  of  the  computation,  one  can 
make  a  rational  decision  in  the  computational  domain.  However,  individual 
familiarity  and  personal  beliefs  often  are  the  dominant  factors  of  to  be 
or  not  to  be  in  one  particular  domain. 


Time-Domain  Calculation  of  Microstrip  Components 
and  the  Curve-Fitting  of  Numerical  Results 


Xiaolei  Zhang  and  Kenneth  K.  Mei 

Department  of  Electrical  Engineering  and  Computer  Sciences 
University  of  California,  Berkeley,  CA  94720 


Abstract 

The  Time-Domain  Finite  Difference  method  com¬ 
bined  with  Fourier  transform  techniques  has  recently 
been  shown  to  be  a  very  effective  approach  in  modeling 
the  dispersive  characteristics  of  microstrip 
coinponents(l-3J.  The  present  research  further  investi¬ 
gate  the  problems  associated  with  the  generation  of 
design  data  over  a  larger  substrate  and  line  parameter 
range  and  the  curve-fitting  of  the  numerical  results  for 
CAD  purpose. 
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A  Super-Absorbing  Boundary  Algorithm  For  Solving  Electromagnetic 
Problems  By  Time-Domain  Finite-Difference  Method 

Jiayuan  Fang*  and  Kenneth  K.  Mei 
Dept,  of  Electrical  Engineering  and  Computer  Sciences 
University  of  California  at  Berkeley,  Berkeley,  CA  94720 

Abstract.  A  new  algorithm  is  presented  to  treat  the  absorbing  boundaries  for  time  domain 
finite  difference  method.  With  this  approach,  the  leading  order  error  of  the  conventional  lo¬ 
cal  absorbing  boundary  conditions  can  be  cancelled  by  a  simple  algorithm,  so  that  the  ab¬ 
sorbing  quality  of  the  boundary  conditions  can  be  greatly  improved. 


TIME  DOMAIN  FINITE  DIFFERENCE  CALCULATIONS 
USING  A  VARIABLE  STEP  SIZE 

Svetlana  Zivanovic*  and  Dr.  Kenneth  K.  Mei 
University  of  California,  Berkeley 


Computational  electromagnetics  (CEM)  is  becoming  increasingly  more  prominent  today, 
solving  problems  which  either  cannot  be  solved  by  analytical  methods,  or  are  too  expensive  to 
investigate  experimentally.  One  such  approach  is  the  finite  difference  method  (EDM) 
developed  by  K.  Yee  (ref.  1)  which  has  been  successfully  used  for  the  calculation  of  elec¬ 
tromagnetic  fields  in  various  structures  (ref.  2). 

However,  one  serious  drawback  of  FDM  is  that  the  computer  memory  required  to  obtain 
sufficiently  accurate  results  approaches  and  often  exceeds  the  available  memory  of  most  com¬ 
puters.  One  way  around  this  is  to  use  bigger  and  bigger  computers  such  as  the  CRAY  II; 
another  is  to  refine  FDM  to  decrease  the  needed  memory  while  keeping  the  error  within  a 
tolerable  limit. 

This  paper  describes  a  new  method  which  uses  a  variable  step  size  (calculation  incre¬ 
ment)  in  the  finite  difference  technique  (VSSFDM).  The  step  size  is  kept  small  around  discon¬ 
tinuities  and  larger  away  from  them,  thereby  using  only  about  1/3  of  the  memory  required  by  a 
straight  application  of  Yee’s  method.  The  results  obtained  by  VSSFDM  are  very  close  to  those 
obtained  by  FDM  (ref.  2)  but  with  a  significant  reduction  in  the  required  computer  memory,  as 
well  as  CPU  time. 


PENETRATION  OF  A  FOCUSED  ELECTROMAGNETIC  PULSE 
INTO  A  BIOLOGICAL  MATERIAL 
AND  ITS  APPLICATION  TO  HYPER  THERMIA 

Hua-Cheng  Chang*  and  Kenneth  K.  Mei 
Department  of  Electrical  Engineering  and  Computer  Sciences 
University  of  California,  Berkeley,  CA  94720 


Abstract  -  A  localized  temperature  distribution  inside  a  biological  material  (human  mus¬ 
cle)  is  obtained  through  the  use  of  a  parabolic  reflector  to  focus  an  electromagnetic  pulse.  With 
the  point-matched  time  domain  finite  element  method,  the  focused  field  and  its  penetration  into 
the  biological  material  are  found;  then,  the  temperature  profile  in  the  material  is  calculated. 
Due  to  the  localization  property,  this  study  can  have  a  significant  contribution  to  the  hyperth¬ 
ermic  therapy.  The  factors  which  influence  the  localization  characteristics  are  also  investigated. 
They  are  the  applied  frequency,  the  associated  electric  parameters  of  the  biological  material, 
and  the  size  of  the  parabolic  reflector. 


Time-Domain  Finite  Difference  Approach 
for  the  Modeling  of  Microstrip  Components 


Xiaolei  Zhang  and  Kenneth  K.  Mei 

Department  of  Electrical  Engineering  and  Computer  Sciences 
University  of  California,  Berkeley,  CA  94720 


Abstract 

The  frequency-dependant  characteristics  of  the 
microstrip  open-end  has  been  analysed  using  several 
full-wave  approaches  [1-5],  The  time-domain  finite 
difference  (  TD-FD  )  method  presented  in  this  paper  is 
another  independant  approach,  which  is  ralatively  new 
in  Us  application  to  obtain  the  frequency  domain  results 
(6].  The  purpose  of  this  paper  is  to  establish  the  validity 
of  the  TD-FD  method  in  modeling  circuit  components 
for  MM1C-CAD  applications. 


ELECTROMAGNETICS  OF  SUPERCONDUCTORS 


K.  K.  Mei,  G.  C.  Liang,  and  T.  Van  Duzer 


Department  of  Electrical  Engineering  and  Computer  Sciences 
and  the  Electronics  Research  Laboratory 
University  of  California,  Berkeley,  CA  94720 


With  the  present  pace  of  research  and  development  activities  in  superconductors,  super¬ 
conductive  boundary  value  problem  is  suddenly  a  “hot”  subject.  To  pursue  this  study,  the 
electromagnetic  behavior  about  a  superconducting  medium  must  be  clearly  understood.  If  a 
superconducting  surface  were  a  lossless  surface  (as  it  is  at  T=0),  it  would  be  electromagneti- 
cally  almost  indistinguishable  from  a  perfecting  conducting  surface  except  that  magnetic  field 
can  penetrate  into  the  superconductor  with  very  small  penetration  depth  (typically  100  nm). 
However,  at  microwave  frequencies  a  superconductor  is  not  a  lossless  material,  and  the  loss  is, 
in  fact,  proportional  to  frequency  square.  Thus,  is  superconductor  just  a  low  loss  conductor? 
What  are  the  boundary  conditions  of  superconductors? 

In  this  paper,  we  will  point  out  that  electromagnetically  speaking  a  superconductor  is 
indeed  distinguishable  from  a  fictitious  perfect  conductor.  A  superconducting  material  is  more 
conveniently  be  treated  as  a  dielectric  material  with  negative  real  part  of  dielectric  constant  as 
far  as  electromagnetics  is  concerned.  By  considering  superconductors  as  a  generalized  dielectric 
material,  some  superconductive  electromagnetic  problems  are  simplified,  since  the  dielectric 
parameter  is  an  integral  part  of  electromagnetic  computation.  It  does  not  present  any  technical 
difficulty  to  existing  computer  programs  if  a  dielectric  constant  passes  from  a  positive  value  to 
a  negative  one. 

In  addition  to  t*  •  conceptional  discussion  of  superconductive  electromagnetics,  this  talk 
will  also  present  a  l  *  domain  computational  method  involving  dispersive  media,  such  as 
superconductors  or  plasma.  The  method  involves  a  time-domain  finite-difference  with  system 
function  expansion,  so  that  no  explicit  convolution  is  required.  It  thus  greatly  reduces  memory 
demand  and  CPU  time. 

As  a  consequence  of  this  electromagnetic  treatment  a  superconductor,  the  existence  of 
a  surface  wave  on  a  superconducting  surface  is  also  predicted.  It  is  a  surface  wave  which  is 
more  tightly  bound  to  the  surface  than  the  conventional  surface  wave. 
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L  INTRODUCTION 

The  use  of  heat  to  cause  elevated  tumor  temperature  at  42°C  or  higher  (hyperthermia) 
has  received  a  great  interest  in  recent  years,  since  in  general  malignant  cells  are  more  sensi¬ 
tive  to  heat  than  normal  cells  in  the  range  of  42-45 3 C  [I].  In  addition,  most  tumors  have 
much  less  blood  perfusion  rate  than  their  surrounding  normal  tissues  have,  meaning  that  they 
may  be  preferentially  heated  [2].  Thus,  obtaining  a  localized  tumor  temperature  has  been  an 
important  task  in  hyperthermic  therapy.  However,  for  non-invasive  treatment,  it  is  still 
difficult  to  produce  well  controlled  elevated  temperature  distributions  (in  the  range  of  42- 
45  °Q  in  the  tumor  without  overheating  normal  tissue.  For  examples,  ultrasound  has  good 
localization  and  deep  penetration  into  human  muscle  but  its  limitation  comes  from  the  high 
absorption  in  bone,  causing  bone  heating.  As  for  electromagnetic  wave,  people  have  been 
using  divergent  waves 'widely  which  often  induce  maximum  heating  on  the  tissue  surface 
with  an  exponential  decrease  with  depth.  Although  low  frequency  does  give  a  good  penetra¬ 
tion,  its  poor  focusing  may  still  hurt  the  normal  tissue  [3].  By  using  a  parabolic  reflector  to 
focus  an  electromagnetic  pulse,  this  paper  presents  a  way  that  can  generate  both  good  focus¬ 
ing  and  satisfactory  penetration,  meaning  that  an  elevated  tumor  temperature  can  be  obtained 
without  overbuming  the  surrounding  normal  tissue. 

The  approaches  in  this  study  contain  the  generation  of  focused  electromagnetic  pulse, 
the  penetration  of  this  pulse  into  a  biological  material,  and  the  calculation  of  the  temperature 
profile  in  the  biological  material.  Pint,  a  plane  wave  of  modulated  Gaussian  pulse  is 
transmitted  towards  a  three  dimensional  parabolic  reflector,  which  causes  a  high  intensity 
field  around  the  focal  point,  as  the  wave  passes.  Next,  this  focusing  field  is  directed  towards 
a  biological  material,  the  time  sequence  of  the  penetration  of  the  fields  into  the  material  can 
be  observed.  The  above  two  tasks  are  completed  by  solving  the  Maxwell’s  equations  with 
the  point-matched  time  domain  finite  element  method  (<i].  Afterwards,  the  temperature 
profile  of  the  biological  material  can  be  calculated  through  the  bio-heat  transfer  equation  [5j. 
Furthermore,  the  localization  properties  of  the  temperature  profile  can  be  adjusted  by  varying 
the  related  parameters. 


This  research  was  supported  by  the  Office  of  Naval  Research  Contract  N00014-86-K-0420. 


-261- 


Time  Domain  Calculation  of  Microstrip  Components 
and  Curve  Fitting  of  Numerical  Results 


Abstract 

The  Time  Domain  Finite  Difference  Method  has  been  used  successfully  to  calculate  the 
Frequency-Dependent  characteristics  of  microstrip  components.  Supercomputers  are  used  to  carry  out 
those  full  wave  solutions.  If  those  data  are  to  be  used  for  MMIC  CAD,  one  needs  to  reduce  them  so 
that  they  can  be  stored  in  P.C.  environment  One  approach  is  to  use  curve  fitting  techniques  such  that 
vast  amount  of  results  can  be  represented  by  appropriate  approximate  functions  with  variable  parame¬ 
ters.  This  approach  appears  to  be  quite  convenient  for  CAD  applications. 


Student  dissertation  resulted  from  this  research: 


Doctoral  Dissertations 

(1)  “Time  Domain  Finite  Difference  Computation  for  Maxwell’s  Equations,”  by 
Fang,  Jiayuan,  December  1988. 

(2)  “Penetration  of  a  Focused  Electromagnetic  Pulse  into  Biological  Material  and  Its 
Application  to  Hyperthermia,”  by  Chang,  Hwa-Cheng,  December  1988. 

Master  Dissertations 

(1)  “On  Radiation  Boundary  Conditions  for  One  Dimensional  Wave  Equations,” 
Master’s  report  by  Fang,  Jiayuan,  May  1987. 

(2)  “Calculation  of  Dispersive  Characteristics  of  Microstrips  Using  Time  Domain 
Finite  Difference  Method,”  Master’s  report  by  Zhang,  Xiao-Lai,  May  1987. 

(3)  “A  Numerical  Analysis  of  a  Shielded  Microstrip  with  a  Finite  Discontinuity,” 
Master’s  report  by  Winnie  P.  Y.  Yu,  August  1988. 

(4)  “Time  Domain  Finite  Difference  Calculations  Using  a  Variable  Step  Size,” 
Master’s  report  by  Zivanovic,  Svetlana,  June  1989. 


